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Abstract 

At variance with the authors' statement [L. Palova, P. Chandra and P. Coleman, Phys. Rev. B 79, 
075101 (2009)], we show that the behavior of the universal scaling amplitude of the gap function in the 
phonon dispersion relation as a function of the dimensionality d, obtained within a self-consistent one- 
loop approach, is consistent with some previous analytical results obtained in the framework of the e- 
expansion in conjunction with the field theoretic renormalization group method [S. Sachdev, Phys. Rev. 
B 55, 142 (1997)] and the exact calculations corresponding to the spherical limit i.e. infinite number N 
of the components of the order parameter [H. Chamati. and N. S. Tonchev, J. Phys. A: Math. Gen. 33, 
873 (2000)]. Furthermore we determine numerically the behavior of the "temporal" Casimir amplitude as a 
function of the dimensionality d between the lower and upper critical dimension and found a maximum at 
d = 2.9144. This is confirmed via an expansion near the upper dimension d = 3. 

PACS numbers: 64.60. an Finite-size systems, 64.60.F- Critical points, 64.70. tg Quantum phase transitions 
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In a recent paper Palova, Chandra and Coleman (PCC)^ studied the quantum paraelectric- 
ferroelectric phase transition in the framework of the self-consistent one-loop approximation ap- 
plied to the familiar quantum ip^ model that plays an important role in the investigations of the 
properties of many quantum systems near their quantum critical points.^ After a suitable nor- 
malization of the parameters of the model this approach is formally equivalent to considering an 
A^-component model in the spherical limit i.e. when the number of the components of the order 
parameter is sent to infinity.^ 

In their equation (33) for the gap function PCC explore the role of the temperature as a boundary 
effect in the "imaginary time". Interpreting the inverse temperature as a finite size in this direction, 
it is possible to map the theory of quantum critical phenomena at low temperatures on the finite 
size scaling theory"^'^ which shaped our current understanding of the modern theory of critical 
phenomena. In this context it is shown in Ref. 1 that at a temperature T above the quantum critical 
point, the gap function A(T) in the phonon dispersion relation scales as 
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The corresponding equation for the scaling amplitude was solved numerically for arbitrary 
dimensions in the range 1 < c? < 3. The bounds on d are imposed by the fact that for d < 1 no 
phase transition can survive, while for > 3 one obtains a mean field critical behavior. Close to the 
upper critical dimension, the presented on FIG. 5 dependence of ad{T 0) on dimensionality 
d exhibits a discrepancy (it is finite) with previous analytical considerations (where it goes to 
zero) as £ ^ 0^ and/or A^ — > oc^, as well as numerical ones for the quantum spherical model.- 
Comparing e-results^ and their numerical prediction PCC suggest that this discrepancy may be 
attributed to the order in which the limits e ^ and A^ ^ oo are evaluated. In this comment 
we demonstrate that at variance with PCC's claims there is no contradiction between previous 
analytical^'^ and numerical^ considerations, and the analysis based on self-consistent one-loop 
approximation presented in Ref. [H 

Our aim is to show that the numerical treatment of Eq. (42) of Ref. jH is not adequate closely 
beneath the upper critical dimension d = 3. For the sake of completeness, we will outline the main 
steps of our computations. According to Eq. (33) of PCC one has: 



where nB{uj) = [c'^/^bT _ i)-i^ = ^n^'^/'^V{d/2), r and 7^ are model constants, and 



\/ q'^ + A^. The parameter = r — Vc measures the distance of the quantum parameter driving 
the transition from its critical value Tc. In the remainder we use kB = 1. Notice that we clearly 
separate the thermal (T > 0) and quantum (T = 0) fluctuations and introduce the cutoff A. 
Further, via the substitution A = aT and q = uT we may write 



T^-M«^_r^ = duu 
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In the low temperature region ^ ^ 1, the cutoff in the first integral can be entirely removed 
neglecting exponentially small corrections. The last integral is convergent in the ultraviolet in 
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dimensions 1 < c? < 3. It may be computed extending the integration over u up to infinity to get 
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where we have introduced the scaling variable 



^'O rp-l/uz 

37c 



with u = (d — 1) ^ the critical exponent measuring the divergence of the correlation length, ^, 
while approaching the quantum critical point and z = 1 the dynamical critical exponent. The 



function is defined via 



u^du 



In particular one has 



yz + v? [exp [yz + v?^ — l] 
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where C,{x) is the Riemann zeta function. 

In the vicinity of the quantum critical point, defined by f^o = and T = 0, the term containing 
in the left hand side may be neglected as it gives only corrections to the leading order of a. 

Then the solution to Eq. ([3]) has the finite temperature scaling form 



ad = fd (x) , 

with fd{x) an universal scaling function. 

For the scaling form ^ to be valid one has to require 



(6) 



,3-d' 



A{d) :-- 



37c 



(47r)'^/2 



r 



d 



0) 



This inequality is an estimation of the low temperature region and the magnitude of 7c where the 
scaling form Q takes place. 

The behavior of the scaling variable above the quantum critical point, obtained numerically 
by equating the right hand side of ^ to zero, is presented in FIG. [B To validate our numerical 
results we choose to perform analytic calculations of ad by considering some particular cases: 
namely d = 2, and in a close vicinity of d = 3 and d = 1. This leads us to the results (see e.g. 
Ref.ii 

'7r(d-l), d-Kl, 
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One sees that the behavior of ad for the cases d ^ 1 and d = 2 seems to be correctly presented 
in FIG. 5 of Ref. 1 . However both numerical (FIG. [B and analytical computations ([8]) show that 
ad vanishes as ^ 3. This result disagrees with the conclusions drawn by PCC based on the 
behavior presented on their FIG. 5. 
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FIG. 1: Dependence of ad = /d(0) on dimensionality d at the quantum critical point i.e. f^o = and 
T 0+. 



In the remainder of this comment we will briefly touch on some aspects of the so called tem- 
poral Casimir effect or Casimir effect in time' considered also some years ago in Ref. in the 
framework of the quantum spherical model. 

The self-consistent one-loop approximation is exact in the spherical limit, i.e. for the theory 
with component order parameter in the limit ^ oo. In this case the free energy can be 
computed using a variational approach based on the Hubbard-Stratonovich decoupling technique- 
of the term (P being A^— component field) in the model of Ref. 1 . Following Ref. 7 we end up 
with an expression for the free energy per particle and per component 



2 sinh ( — C(j, 
2T 
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The equation minimizing the free energy ^ is identical to the self consistency equation ([T]) up to 
the substitution 7c 37c. Notice that expression ^ for the free energy can be split into a "pure 
quantum" free energy and a "finite temperature" contribution as: 

^^(^' -Ji" {-4 - '■)^ n / (g^"' + ^ / S (1 - -"•'^) } • <«» 

Applying the assumption of Ref. 10 for the singular part of the free energy for classical sys- 
tems within the theory of finite size scaling to the case of quantum critical phenomena at finite 
temperature, the singular part of the free energy should scale like 

where the function ^ is an universal scaling function, whose value at the quantum critical point 
is equivalent to the Casimir amplitude^ for temperature driven phase transitions in films. Such 
an idea, considering the inverse temperature as an additional dimension, has been developed for 
quantum critical points in Ref. 9. Notice however that for quantum systems this quantity may be 
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measured experimentally since it is related to the amplitude of the specific heat of the system at 
finite temperature. 

For the model under consideration in the interval of interest i.e. dimensions 1 < c? < 3 and 
in the vicinity of the quantum critical point, the singular part of the free energy takes the scaling 
form 

^fs-(r) = ri+^^7d(x), (12) 

where 

9.{>^) = \-o.l - i(4.)-('^+^)/^r - n-^'^'^/'T h,^M) (13) 

is an universal scaling function with ad the solution Q to Eq. ([3]). It is worth mentioning that the 
scaling form (fT2)) is in agreement with the scaling ansatz (fTT)) . 

The behavior of the scaling function (fT3l) for d = 2 is well known in the literature.-2iii Here we 
wish to check its dependence upon the dimensionality. We are primarily interested in the behavior 
of the amplitude gdiO) of the free energy at the quantum critical point, as this is tightly related 
to the "Casimir effect in time".- For the particular case d = 2, it can be computed analytically 
resulting in^ 

^2(0) = (14) 

For arbitrary d the behavior of (7^(0) can be obtained by numerical means. This is graphed in FIG. 
m It is found that the scaling function has a maximum at d = 2.9144. We check our numerical 
results in the vicinity of d = 3, using an e = 3 — d expansion, taking into account the small e 
behavior of as-^, to get 
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where 7 = 0.5772 is the Euler-Mascheroni constant. Indeed, expression (fT5l) exhibits a maximum 
atd = 2.9818. This maximum is shifted compared to the one obtained by numerical means due to 
the used approximation for e. 

At the borderlines d ^ 1^ and c? ^ 3"*", the values of gd{0) coincide with their counterparts of 
the Gaussian theory corresponding to = i.e. 

^Gaus.^Q) = _^-(d+l)/2r (^^^ C(rf+ 1) 

This can be seen in the inset of FIG. |2l where we present the behavior of 

P = ^7d(0)/^?f^^^-(0). (16) 

which is related to the Zamolodchikov C-f unction extended to arbitrary dimensions and nonzero 
temperature (see Ref. S and references therein). Here, we would like to point out the similarity 
between the behavior of p in FIG. [Hand the one obtained in Ref. llJfor the characteristic parameter 
of the corresponding conformal field theory in dimensions 2 < c? < 4. 

Let us note, before closing this comment, that the gap equation [Eq. (|3])] with l.h.s. equals zero 



is equivalent to the spherical constraint imposed on the quantum spherical model, see Refs. I8ll9lll3 
where the finite-temperature scaling was studied in great details. 
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FIG. 2: The behavior of the Casimir amplitude gd{0) from Eq. (fT3l) . In the inset we graph the ratio p 
defined in Eq. (fT6l ). 

The consideration outlined in this comment remains valid also for systems with film geometry 
under periodic boundary conditions with temperature driven phase transition. This is due to the 
fact that the thickness of films is in some sense equivalent to the inverse temperature in a quantum 
system. This is a facet of the property known as temperature inversion symmetry discussed in the 
literature that leads to explicit conversion from Casimir force to Planck's law of radiation.-^^ Very 
recently, in Refs. 16 and 17, using the e expansion and/or the limit ^ oo in the framework of 
the classical 0{N) symmetric Lp^ model with film geometry one obtains results for the amplitude 
of the correlation length and quantities related to the Casimir effect in a close vicinity of the upper 
critical dimension that can be conversed to the field of quantum paraelectric-ferroelectric phase 
transitions in particular and to quantum critical phenomena in general. 
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